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Abstract 

Putting dynamics into random matrix models leads to finitely 
many nonintersecting Brownian motions on the real line for the eigen- 
values, as was discovered by Dyson. Applying scaling limits to the 
random matrix models, combined with Dyson's dynamics, then leads 
to interesting, infinite-dimensional diffusions for the eigenvalues. This 
paper studies the relationship between two of the models, namely the 
Airy and Pearcey processes and more precisely shows how to approxi- 
mate the multi-time statistics for the Pearcey process by the one of the 
Airy process with the help of a PDE governing the gap probabilities 
for the Pearcey process. 

1 Introduction 

Putting dynamics into random matrix models leads to finitely many nonin- 
tersecting Brownian motions on K for the eigenvalues, as was discovered by 
Dyson [13]. Applying scaling limits to the random matrix models, combined 
with Dyson's dynamics, then leads to interesting, infinite-dimensional diffu- 
sions for the eigenvalues. This paper studies the relationship between two of 
the models, namely the Airy and Pearcey processes and more precisely shows 
how to approximate the Pearcey process by the Airy process with the help 
of a PDE [2] governing the gap probabilities for the Pearcey process. The 
Airy process was introduced by Prahofer-Spohn |22j and further developed 
by K. Johansson [T71 ITS] . A simple non-linear 3rd order PDE for the tran- 
sition probabilities for this process was found in [3] ; see also [231 [21] . The 
Pearcey process was introduced in [231 EI] in the context of non-intersecting 
Brownian motions and plane partitions, also based on prior work on matrix 
models with external source [2Ql HDl ESI EZl HH H31 El [3 El E] - 
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Consider n nonintersecting Brownian particles on the real line R, 



— OO < X\{t) < ... < x n {t) < oo 
with (local) Brownian transition probability given by 

p(t;x,y) := 



nt 

all starting from the origin x = at time t = and such that | particles 
are forced to ±y^f at £ = 1. For very large n, the average mean density of 
particles has its support, for t < |, on one interval centered about x = 0, 
and for | < £ < 1 on two intervals symmetrically located about x = 0. The 
end points of the interval(s) of support describe a heart-shaped region in 
(x, i)-space, with a cusp at (xo,t ) = (0,1/2). The Pearcey process V(t) 
is defined (see Figure 1) as the motion of these nonintersecting Brownian 
motions for large n, about (a;o^o) = (0, 1/2) (i.e., near the cusp), with space 
microscopically rescaled by a factor of n -1 / 4 and time rescaled by a factor 
n -1 / 2 , in tune with the Brownian motion rescaling. A partial differential 
equation for the Pearcey process was found by Adler-van Moerbeke jl] and 
in [2] a much better version was obtained, namely a simple third order non- 
linear PDE for the transition probabilities; it was obtained by a scaling limit 
on a PDE for non-intersecting Brownian motions with target points. This 
PDE is related to the Boussinesq equation and its hierarchy; this is part of 
a general result on integrable kernels, as explained in the forthcoming paper 

®- 

Near the boundary of the heart-shaped region of Figure 1, but away from 
the cusp, the local fluctuations behave as the so-called Airy process, which 
describe the non-intersecting brownian motions with space stretched by the 
customary GUE edge rescaling n 1 / 6 and time rescaled by the factor n 1 / 3 , 
again in tune with the Brownian motion space-time rescaling. 

This paper shows how the Pearcey process statistics tends to the Airy 
process statistics when one is moving out of the cusp x = J^(3(£ — t )) 3 / 2 
very near the boundary, that is at a distance of (3T) 1 / 6 for r very large, with 
t being the Pearcey time. To be precise, in the two-time case, the times Ti 
must be sufficiently near -in a very precise way- for the limit to hold. The 
main result of the paper can be summarized as follows: 
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Theorem 1.1 Given finite parameters t\ < t 2 , let both T\, t 2 — > oo, such 
that T 2 — T\ — > oo behaves in the following precise way: 




T<1 ~~ Tl (Q M/3 , ^2 — ^1 2tit 2 ^/ -5/3\ /-, -i \ 

^^) = (3ri) + (3^ + ^T + ^ )• (L1) 

The parameters ti and t 2 provide the Airy times in the following approxima- 
tion of the Airy process by the Pearcey process: 



( 2 \ / \ (L2) 

= p n wo n (-^) = }j f 1 + 0(773) j • 

The same estimate holds as well for the one-time case. A similar (but 
different) result was then obtained in [6J for the one-time case in the situation 
when one is moving out of the cusp following both the two branches of the 
cusp simultaneously. 

The proof of this Theorem proceeds in two steps. At first, Proposition 



2.1[ stated later, shows that, using the scaling above, the Pearcey kernel 
tends to the Airy kernel. In a second step, we show, using the PDE for the 
Pearcey process [2] and Proposition 2.1 the result of Theorem |1.1| 



It is a fact that both the Airy and Pearcey processes are determinantal 
processes, for which the multi-time gap probabilities is given by the matrix 
Fredholm determinant of the matrix kernel, which will be described here. 

The Airy process is a determinantal process, for which the multi-time 
gap probabilities are given by the matrix Fredholm determinant for t\ < 
...<t e , 

p ^Q(^) n E * = )) = det ( J - ^ K A^) 1<lJ<e 

of the matrix kernel in x = (xi, . . . ,xe) and y = (yi, ■ ■ ■ ,yt), denoted as 
Xe is the indicator function for the set E 
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follows: 



I KfifrfcuyxWdxidyi ■■■ K£,t t ( x uye)\ / dx 1 dy e \ (1.4) 

\ K h,H ( x i, Vi)Vdxedyi ■ ■ ■ Kf tM (x e , y t ) yjdx t dyz ) 

where, for arbitrary ti and tj, the extended Airy kernel K A t .(x,y) is given 
by (see Johansson |TTl ITS] ) 

K& tj ( x , V) ■= K£ t . (x, y) - 1(U < t j )p A {t j - t h x, y), (1.5) 

where 

POO 

K A t .{x,y):= / e- x ^A(x + X)A(y + X)d\ 
Jo 

Iff e -v s /3+yv i 

= {2mf J T> V Jr < U e- 3 /3+^ (v + tj ) - (u + U) 



p A (t,x,y) := . l " ^ K y+X \ for t > 0. 

V47Tt 

The contour k G T< consists of two rays emanating from the origin with 
angles 9\ and 8[ with the positive real axis, and the contour v G T> also 
consists of two rays with angles 82 and #2 with the negative real axis, as 
indicated in Figure 2. As is well known, one may choose 8\, 82, 8[, 8' 2 = ir/3. 
In particular, for tj = tj, this is the customary Airy kernel 



POO 

K A ti := K A (x, y) := / d\A(x + \)A(y + A) 

Jo 



A(x)A'(y) - A{y)A\x) 
x — y 

(1.6) 



4 



n/2 # of particles n/2 




x 

Figure 1: The Pearcey process. 




3tt/8 < t,t' < 5?r/8 
Pearcey contour Airy contour 

Figure 2: Integration paths for the kernels. 
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The Pearcey process is determinantal as well, for which the multi-time 
gap probability is given by the Fredholm determinant for n < . . . < Te, 



for the Pearcey matrix kernel in £ = (£1, ...,&) and fj = (r]\, 
denoted as follows: 



Tl,...,T£ 



;i-7) 



\K% ! T 1 (&>Vi)\/d&d'ni ... K% tTt {€ e ,ri e )y/d£ e dri e J 
where for arbitrary Tj and Tj, (see Tracy- Widom [25J) 



< >T . (£, /y) = if' (£, 77) - 1(ts < ^(r, - r,; £, ,7) 



-7> 



V, 



with 



1 



K Utn):=-^j x du 



dV 
Y V-U 



1 T i v 

e 4^2 



-Vr) 



e 4 



-e 2t , tor r > 0, 



;i.9) 



\/27rr 

where X and Y are the following contours: X = X\ U X2 consists of four 
rays emanating from the origin with angles a\, a[ with the positive real axis 
and <J2,<j' 2 with the negative real axis, as given in Figure 2. The contour 
Y consists of two rays emanating from the origin with angles r, t' with the 
negative real axis; it is customary to pick a\ = 02 = o' x = a' 2 = tt/A and 

T = T' = 7T/2. 

In [U |2], it was shown that, given intervals 



the log of the probability 



Q = Q( n , . . . , Tf, E x , . . . , E t ) := logP f]{V(r t ) n E t = 0} 



vi=l 



satisfies a non-linear PDE, which we describe here. Given the times and the 
intervals above, one defines the following operators: 

^ := E|:' ^ := Ei 5 -=E^ 

i 1 k "Sfe i 

*-E*£. E « : =EEf?i ("») 

Then Q satisfies the Pearcey partial differential equation in its arguments 
and the boundary of the intervals: 

2d 3 T Q+\(2e T +e E -2)d 2 E Q-(J2nd E Jd T d E Q+{d T d E Q,d 2 E Q} g =0. (1.11) 



2 From the Pearcey to the Airy kernel 

Define the rational functions 

1 t x - t 2 4 U o 

f ; u) := - - -— g + — + -tx + -t x 

4(3m) j {ouy 3u 3 6 ^ 1) 

h(x,t;u) := — (x + 6t 2 ), 



and the diagonal matrix 
S : = 



e ^(x,t 1 ;z 4 )-h(x,t 1 ;z 4 ) q 



We now state: 



Proposition 2.1 Given a parameter z — > 0, define two times T\ and r 2 
blowing up like z~ % and depending on two parameters t\ and t<i, 

r l = ^(l + 6^ 4 ) + 0(^ 1 °). (2.2) 
Given large ^ and rjj, define new space variables Xi and ijj, using 7$ above, 

= hln?' 2 - (Sn)^ 

27 (2.3) 



With this space-time reseating, the following asymptotic expansion holds for 
the Pearcey kernel in powers of z A , with polynomial coefficients in t\ + t 2 , 



S K^^,ff)y/d^dff S- 1 = (1 + (t 1 +t 2 )0 1 (/)) Kf ut2 {x,^dxdy + 2 (z s ), 

(2-4) 

where 0\ refers to a differential operator in d/dx,d/dy, with polynomial 
coefficients in x,y,t\ ± t 2 and (t\ — ^2) 1 and varying with the four matrix 
entries of the matrix kernel. 

On general grounds, due to the Fredholm determinant formula, this would 
give Theorem but with a much poorer estimate in (1.2), namely instead 
of 0(ti ), the estimate would be 0(t(( 2 ^ 3 ). The existence of the PDE for 
the Fredholm determinant of the Pearcey process enables us to bootstrap the 
0(t x 2 ^ 3 ) estimate to 0{r t 4 ^ 3 ), without tears! 

Also note that conjugating a kernel does not change its Fredholm determi- 
nant. Before proving Proposition 2.1, we shall need the following identities: 

Lemma 2.2 Introducing the polynomial 

1 3 

s; u) := -u 4 + ^uj 2 s 2 - As(xuo - w 3 ), (2.5) 



the Airy kernel (1.6) and the (double) integral part of the extended Airy kernel 



(1.5), at t\ — s, t 2 = —s, satisfy the following differential equations, 
8, -d x ) - V(y, s, d y ) + 4s) K A (x, y) 



-{x - y){x + y + 6s 2 )K A (x,y), 



(2.6) 



and 



U(x, s, -d x ) - V(y, -s, d y ) - \s^{d x - d y )j K A _ s (x, y) 



\{x - y){x + y + 6s 2 )K A s (x, y). 



(2.7) 



Proof: The operator on the left hand side of (2.6) reads 
^{(x,s,-d x )-^{(y,s,d y )+As 

= 4,(1 + yd y - dl + xd x - dl) + \s\dl - d 2 y ) + \{d' x - d A ). m 



Using the first representation (1.6) of the Airy kernel, one checks using the 
differential equation for the Airy kernel, A"(x) = xA(x) and thus A"'(x) = 
xA'(x) + A(x) and A^ iv \x) = 2A'{ x) + x 2 A(x), and using differentiation by 
parts to establish the last equality, 

(( ydy -d 3 y ) + (xd x -dl))K A (x,y) 

= - J dz (^z-^ + 2J A(x + z)A(y + z) = -K A (x,y). 



In order to take care of the other pieces in (2.8), one uses the second repre- 
sentation (1.6) of the kernel K A (x,y), yielding 

(dl-d*)K A (x,y) = (x-y)K A (x,y) 
{dt-tf)K A (x,y) = (x 2 -y 2 )K A (x,y). 



This establishes the first identity (2.6) of Lemma 2.2 The operator on the 
left hand side of the second identity (2.7) reads: 

3 d 

s, -d x ) - $?(y, -s, d y ) - 2 S ^(^' _ d v> 



3 d 



J%-9S) + + - di - yd y + di)) - - s -(d x -d y ), 

of which we will evaluate all the different terms acting on the integral. Notice 
at first that, since the expression under the differentiation vanishes at and 
oo, one has 



roc a 

= y dz—(ze- 2sz K A (x + z,y + z)) 

POO POO 

= / dz e~ 2sz K A (x + z,y + z) - / zdz 
Jo Jo 

zdz e~ 2sz K A (x + z,y + z), 



-2sz 



A{x + z)A{y + z) (2.9) 



Again, using the second representation (1.6) of the kernel K A (x,y), one 
checks 

POO 

(d x -d y )K A _ s (x,y) = -(x-y) / dz e~ 2sz K A {x + z, y + z) 

Jo 

(dl - d 2 )K A _ a (x,y) = (x- y)K A _ s (x,y), 



and also 

d f°° 
s — (d x -d y )K A _ s (x,y) = 2s(x - y) j zdz e~ 2sz K A (x + z,y + z 

and, using the differential equation (xd x — d x )A(x) = — A(x) and (2.9), 
-2s {{xd x - SQ - (yd y - $)) Kf_ s {x,y) 

POO 

= -2s(x-y) zdz e~ 2sz K A (x + z,y + z) 



Using AW(x) = 2A'( x)+x 2 A(x) and the Darboux-Christoffel representation 
( |1.6[ ) of the Airy kernel, one checks 

(%-ai)kf_,(x, y ) 

e~ 2sz (A^ v \x + z)A{y + z) - A{x + z)A (iv \y + z))dz 

POO 

= -2{x-y) / dz e~ 2sz K A (x + z,y + z) + {x 2 - y 2 )K A s (x,y) 
Jo 

POO 

+2(x-y) zdz e~ 2sz A(x + z)A(y + z). 
Jo 

Adding all these different pieces and using the expression for 

POO 

—2s(x — y) / zdz e~ 2sz K A (x + z,y + z), 
Jo 



given by (2.9), leads to the statement of Lemma 2.2 



For the sake of notational convenience in the proof below, set 

^((blilvW Kf 2 (^i,r} 2 )Vd^drf 2 

^21(6^1)^2^1 #22 (6,^2) Vd&chfr 



K^mdidfj 



(2.10) 



and similarly for the Airy kernel K^ t2 ; the and K A refer, as before, to 
the (double) integral part. 



Proof of Proposition 2. 1 : Notice that acting with d x and d y on the kernels 

K A (x, y) and K A (x,y), amounts to multiplication of the integrand of the 

22 21 
kernels with —u and v respectively; i.e., u -H- —d x and v -H- d y . 
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Given the (small) parameters e R, consider the following (z, j)-dependent 
change of integration variables ([/, V) i— > (u, v) in the four Pearcey kernels 
Kfj in (2.10), namely 



U= -^(l + 3^ 4 )(l + 3t^ 4 ), V= -^(l + 3^ 4 )(l + 3t^ 4 ), (2.11) 
together with the changes of variables (£, 77, ts, Tj) (x,y,ti,tj), in accor- 



dance with (2.2) and (2.3) 
1 



7* = —11 + 6^ 4 ) + 0(z lu ), = —(1 + 6t^ 4 ) + 0(2 

oZ oZ 

Z = |(3r,) 3/2 - (Sr,) 1 / 6 *, r, = |(3r,) 3 / 2 - (St^V 



10 \ 



(2.12) 



To be precise, for K]? k , one sets in the transformations above i = j = k, 
for K^ 2 and K^-, one sets i — 1, j = 2 and z = 2, j = 1 respectively. 



Then, remembering the expressions $ and defined in (2.1) and (2.5), the 
following estimate holds for small z: 



e -*(v,*i;**) — ^ e -« 4 *(2/,*i^) e -\+yv 



l+^0(z 8 )+^-0(z 8 )) 



e 3 



e -z 4 9(x,ti;-d x ) ^_y? +xu 



l+UOiz^+tjOiz 8 )); 

(2.13) 

replacing in the latter expression u and f by differentiations d x and <9<, has 
the advantage that upon doubly integrating in u and v , the fraction of ex- 
ponentials in front can be taken out of the integration. Moreover, setting 
t\ = t + s and t 2 = t — s, one also checkg^] : 



e -$(a;,t i ^ 4 ) 



V 



v 



p A (—2s, x, y) a/ dxdy ( 1 + — [(x— y)(x+y + 6s 2 ) H — r(x, y, s)] + 0(z 8 ) 



(2.14) 



5 The precise expression for r(x, y, s) := (x — y) 2 — 8s 2 (x + y — 2s 2 ) will be irrelevant in 
the sequel. 
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The multiplication by the quotient of the exponentials on the left hand side 
of the expressions above will amount to a conjugation of the kernel, which 
will not change the Fredholm determinant. Setting t * = t + s and tj = t — s, 
with t — 0, one finds for the non-exponential part in the kernel 



dUdV 



v-u 

yjdxdy dudv(l + |(t; 



v + tj — u — ti + 3z 4 (vtj 



Uti 



+ 0{z i 



V 'dxdy dudv 
v—u 



l±4s2 4 + 0(z 8 )) for 



V 'dxdy dudv I -y | 3z 4 s(u+v) 



+ 0(z 8 



for 



1. 
2- 

i 
i 



i = l 

/ = 2 

1, 3 

2, j 



2 
1. 



(2.15) 



Along the same vein as the remark in the beginning of the proof of 
this Theorem, one notices that multiplication of the integra nd of the kernel 
K A (x, y) with the fraction, appearing in the last formula of (2.15), amounts 

21 I 

to an appropriate differentiation, to wit: 



, ?>z A s(u + v) 3 d 

± < — > -s — (d v 

v-ut2s 2 ds K y 



(2.16) 



So we have, using (I2.13L 12.141), (|2.15|), and the differential equation for K A 



of Lemma 2.2 



i ;2 4 ) 



K' A (x, y^\fdxdy 



t=o 



(±4s + ^((x, ±s, -^.) - *((</, ±s, 9„)) i^(x, y) + 0(z 8 ) 

(2.17) 



x 



y)(z + y + Qs 2 )K A (x, y) V 'dxdy + 0(; 



12 



12 
21 



2.2 



the upper (lower)-indices correspond to the uppe r (low er)-signs. Using (2.16), 
and the differential equation for K A of Lemma 



-$(x,t i ;z 4 ) 



ttzi\ K v 12 {£, rf) ^d^drj - Kf (x, y) ^/ dxdy 



t=o 



3 d 



+ 0(z 8 



(x - y) (x + y + 6s 2 )K A (x, y) \fd~xdy 



t=o 



Also, using (2.14), 



-*(y^2^ 4 ) 



.q,(^ijP( T 2 ~ n;£,v)V d £ dr l -p A (~'2s,x,y)^dxdy 



i=0 



Z 

T 



(X 



- y)(x + y + 6s 2 )p A (-2s, x, y) a/ dxdy + 0(z 8 ). 



(2.18) 



(2.19) 



Since h(y, ti, z 4 ) — h(x, tj] z 4 ) 
at t = 0, and thus 

-h(y,ti;z 4 ) 



^-(x—y)(x+y+Qs 2 ) for arbitrary Uj = t±s 



1 + — (x -y){x + y + 6s 2 ) + 0(z s ) 



(2.20) 



e -h(x,tj;z 4 ) 1 4 

Then the following approximations follow upon combining the three estimates 
above d2l7h, (12181), d2ll and using d^20h: 



-*(j/,ti;« 4 )+fe(l/,ti; J i 4 ) 



-*(a:,t i ;2 4 )+/i(:r,t i ;z 4 ) " 11 

e 2 2 22 



- K A (x,y)^dxdy = 0(z 8 



t=o 



-^{y,t 2] z 4 )+h(y,t 2 ;z i ) 



-$(x,t i ;z 4 )+h(x,t i ;z 4 ) 12 
g 2 2 21 

-$(2/,t2;2 4 )+My,t2;z 4 ) 



-K A (x,y)^dxly = 0(z 8 ) 



-*(x,ti;2 4 )+fc(x,ti;z 4 ) 



t=0 



t=0 



p A (-2s,x,y)^dxdy = 0{z 8 ). 



The reader is reminded that this estimate sofar is done at t — 0. It then 
follows for t ^ from the estimates (pill, ((2121), ((2131), ((211, plB], 
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(2.16) that the right hand side of (2.4) is an asymptotic series in z A with 
polynomial coefficients in t\-\-t 2 = 2t, proving the claim about 0\. 

So far an important point was omitted, namely to analyze how the Pearcey 
contour turns in the limit into an Airy contour; a detailed description of the 
possible contours was given in Figure 2. The Pearcey-rays X±, with an- 
gles <J\-,o'\ G (7r/8,37r/8) can be deformed into two acceptable rays 9i,6[ G 
(7r/6, 7r/2) for the T<-Airy contour, since the two intervals have a non-empty 
intersection. Also, since (37r/8, 57r/8)n(7r/6, 7r/2) ^ 0, the Pearcey Y-contour 
can be deformed into an acceptable v G r>-Airy contour. Again, since the 
admissible interval (7r/8,37r/8) for cr 2 and a' 2 in the X 2 -Pearcey contour con- 
tains 7r/6, one ends up in the limit integrating the function e u / 3 along a 
contour of the form r>, which we may choose to have an angle of 7r/6 — 5 
with the negative real axis, for small arbitrary 5 > 0. In the sector r>, 
with < 9 2 = 9' 2 < 7r/6 — 5, the function e u / 3 decays exponentially fast. 
Therefore, the contribution, due to the it-integration of e u / 3 x (lower order 
terms) over r>, having an angle of 7r/6 — 5 with the negative real axis, van- 
ishes by applying Cauchy's Theorem in that sector. This ends the proof of 



Proposition 2.1 



3 From the Pearcey to the Airy statistics 



This section concerns itself with proving Theorem |1.1| 

Proof of Theorem \l.l\ For any intervals E% and E2, the function 

Q(t 1 ,t 2 ;E 1 ,E 2 ) = logP (n(7>(T f ) n E { = 0) ) (3.1) 



satisfies the Pearcey PDE ( 1.11 ). Reparametrizing, without loss of generality, 



Tl = T + a, r 2 = T-a, Ei = (£ + r] + fi,£ + ri-n), E 2 = (£-il + v,£-r]-v), 

(3.2) 

leads to a manageable PDE for the function 

F(r, a; £, r], //, v) := Q(t u t 2 ; E u E 2 ), (3.3) 
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namely the PDE: 



d 3 F 1 
+ 4 



2{ a —-T — )+i—+ — + — + is— -2 s °' 2 ''' 
° da T dr <9£ ^ drj d^i V dv 



d 2 F f d 2 F d 2 F 



(3.4) 



0. 



To go from the Pearcey PDE ( 1.11 ) to the PDE (3.4), one notices that the op- 
erators (1.10), appearing in the Pearcey PDE (1.11), have simple expressions 
term; 



in terms of the variables r, a, £, i], /i, z/, 
dF „ „ dF 



dr 



d E Q 



^ />, d d d 3. 



The change of variables, considered in (|2.2[ ) and (2.3), combined with a 
linear change of variables, in parallel with (3.2), 



Ti 



E: 



1 

3? 
2 
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1 + Qt t z 4 ) with 
(3nf 2 - (Zr^Ei with 



ti=t + S 

t 2 =t - s 

Ei = (x + y + u, x + y — u) 
E 2 = (x — y + v , x — y — v ) 



(3.5) 



yields a z-dependent invertible map, 

T: (t,s,x,y,u,v)i-> (t,<t,€,7i,h,v) 



(3.6) 



and thus the function F in (3.3) leads to a new z-dependent function G: 
F(t, a; f , ij, fi, u) = F(T(t, s, x, y, u, v)) =: G(t, s, x, y, u, v), 



of which we compute, in principle, the series in z. From Proposition |2.1[ it 
follows that the z 4 -term in the asymptotic expansion in (2.4) vanishes when 
t — 0; so, omitting ^dxdy, one has 



SK^S- 1 = Kf ut2 + z 4 Kx + Z ^ K ^ with K i = m > 
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for some kernel H, with the Kj polynomial in t. Then defining 



(3.7) 



one finds: 



G(t,s,x,y,u,v) 

2 

= logP(p|{P(r i )n£ l = 0} 



i=l 



logdet(/-< )T2 ) EiXj32 
logdet(I-K£j. - 



{3A 



z a TtL 1 -z 8 Tt(L 2 + -L 2 1 ) 



G (s; E 1} E 2 ) - G^t, s; E u E 2 )z 4 - G 2 (t, s; E u E 2 )z s + 0(z 12 ) 



where 



G (s; E U E 2 ) = log det(J-K^ 2 



«1 X E 2 ' 



(3.9) 



Gi{t, s; E u E 2 ) = t Tr((7 - K^J^H) 



XB 2 



Note Go(s;Ei,E 2 ) is t-independent, since the Airy process is stationary. 
Then setting 

F(r, a; £, rj, v) = G(T _1 (r, a; f , r], //, z/)) = G(t, s, x, y, u, u) 



into the PDE (3.4) yields, using differentiation by parts, a new PDE for 
G(t, s, x, y, u, v). Upon setting the expansion (3.8) of G for small z, 

G(t, s, x, y, u, v) = G (s; E u E 2 ) - G 1 (t, s; E x , E 2 )z 4 + 0(z 8 ), 

into this new PDE leads to 

d 3 G d 3 d 



z z 2s 



dsdx 2 dtdx" 2 



+ 0(z b 



0. 



(3.10) 



with G\ being polynomial in t. From the term in the PDE (3.4), it 
would seem like the leading term would be of order z~ 6 ; in fact there are 
two consecutive cancellations, so that the first non-trivial term has order 
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z , thus leading to a simple PDE connecting Gi with Gq. But since Go is 
t-independent, equation (3.10) can also be written as 

d d 2 { d \ 

m* G '- 2t ^ G ° = ° (311) 



From (3.9), one finds 



0_ 

ds 



G 1 - 2ts^-G Q = t ( Tr((J - K^J^H)^ - 2s — G 



' E 1 xE 2 



d_ 
ds 



^2t l ai(s;x,y,u,v), 



which substituted back in (|3.11|), leads to the PDE's for the a^, namely 

ai(s;x,y,u,v) = 0, (3.12) 



( d 



- 



implying 



\dx 



CLi(s; x, y, u, v) = xbi(s; y, u, v) + q(s; y, u, v). 



(3.13) 



Letting the intervals E\ and E2 go to 00, while keeping their relative position 
2y fixed and wid ths — 2u and —2v fixed as well, is achi eved by letting x — > 00, 
as follows from (3.2). Remember t±,t2, Ex, E2 from (3.5). But in the limit 
x y 00 , the expressions G , namely 

G (s;E u E 2 ) = \ogdet(I-Kl t2 )- ix&2 

tends to exponentially fast using the exponential decay of the four compo- 
nents of the matrix Airy kernel (1.5); note that the term p A (t,x,y) tends to 
as well, when x and y tend to 00, due to the presence of e~ t( - x+y ^ 2 . Next, 
we sketch the proof that 



G 1 (t, s; Ei, E 2 ) = t Tr((J - K^J^H) ^ 



XE 2 



tends to exponentially fast, when x — > 00. Indeed, using the identity (R is 
the resolvent of the Airy kernel t ) 



(J - Kf^M = H + Kf ut2 (I - K^J^H = H + Kf ut2 (I + R)H 



-A 



A 



-A 
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one computes 

G 1 = 'ft((J-K^ >lsl )- 1 K 1 ) 
= tTr((I-K? ut2 y l M) &i> 

= tT ^ 1 ^ 2 + tT ^tfi) &lX ^ with H := (/ + R)H, 



E 1 xE 2 
L xE 2 

-A 



where 



2 r 

TrH ^ 2 =E/~ Hu(u,u)du 

1=1 JE i 

2 

Tr(Kf 1>t2 M)- ix - 2 = fi . K£ ti (u,v)H u (v,u)dudv 



i=l 



+ // K£ t2 (u,v)H 2 i(v,u)dudv 

J J EixE 2 

+ // K^ tl (u,v)Hi 2 (v,u)dudv. 

J J E 2 xE± 

Each of these integrals tend to because each of them contains the Airy 
kernel and also because H is obtained by acting with differential operators 
on the Airy kernel, as explained in section [2] This ends the proof that 
G\(t, s; Ei, E 2 ) —> 0, when the Ei tend to oo. This fact together with the 
form (3.13) of the dj imply b{(s; y, u, v) = q(s; y, u, v) = for i > 1, and thus 
di = for i > 1, implying 

d 

Gi = 2ts— Gq. 

OS 

Summarizing, this implies that for Ei = ^(3rj) 3//2 — [?>Ti) l ^Ei, substituting 
Gi = 2ts(dG /ds) into §M, 



\i=l 



G (s;Ei,E 2 )-2zHs—^ + O(z 8 ) 

OS 

G (s-2tsz 4 ;Ei,E 2 ) + O(z 8 ) 
logP (q - Uz 4 )) n {-Ei) = 0} ) + 0(z* 
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in view of the change of variables given in (3.5), one has that s — 2tsz 4 = 
§(£i(l — tiz A ) — £2(1 — t 2 z 4 )). Then, substituting £j = + u,iz 4 ) into 
U{1 - Uz A ) and 75 = (1 + 6t iZ 4 )/(3z 6 ) + 0(z 10 ), as in Q, yields 

1 2m- 

ti(l-t^ 4 ) = Ui -2uh s -u 4 z 12 and Ti = — + ^+2m^ 2 +0( z 10 ) for i = 1,2. 

3z b 

Then eliminating z between the two expressions above for t\ and T2, by first 
expressing 2 as a series in t±, yields (1.1) and (1.2), with U replaced by Uj. 



This ends the proof of Theorem 1.1 
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